Nonlinear phase-dynamics in a driven Bosonic Josephson junction 
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We study the collective dynamics of a driven two mode Bose-Hubbard model in the Josephson 
interaction regime. The classical phase-space is mixed, with chaotic and regular components, that 
determine the dynamical nature of the fringe-visibility. For weak off-resonant drive, where the 
chaotic component is small, the many-body dynamics corresponds to that of a Kapitza pendulum, 
with the relative-phase ip between the condensates playing the role of the pendulum angle. Using a 
master equation approach we show that the modulation of the inter-site potential barrier stabilizes 
the tp — n 'inverted pendulum' coherent state, and protects the fringe visibility. 



The Josephson effect [l[ is an unambiguous demonstra- 
tion of macroscopic phase coherence between two cou- 
pled Bose ensembles. The experimental realization of 
BEC Josephson junctions [2™ has led to observations of 
Josephson oscillations 0, 0, US] as well as macroscopic 
self-trapping 0, B ancl the equivalents of the ac and dc 
Josephson effect [1, [l(| , present also in the superconduc- 
tor [12| or the superfiuid Helium Lj, [lj] realizations. 
Beyond these mean-field effects, BEC Josephson systems 
allow for the observation of strong correlation phenom- 
ena, such as the collapse and revival of the relative phase 
between the two condensates 15l4l7| which was observed 
with astounding precision in an optical lattice in Refs. 
fl8| . in a double-BEC system in Ref. and in a ID 

spinor BEC in Ref. 

The bosonic Josephson junction is often described by 
a two-mode Bose-Hubbard Hamiltonian (BHH) [2l|, [22[ , 
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Here K, E, and U are coupling, bias, and interaction en- 
ergies. The SU(2) generarators L x = (a\a2 + a\ai)/2, 
L y = (a\a.2 — a\a\) / (2i) , and L z = (n\ — ri2)/2, are de- 
fined in terms of the boson on-site annihilation and cre- 
ation operators hi, fit, with the conserved total particle 
number ri\ + 112 = N = 21. Three distinct interaction 
regimes are obtained, depending on the characteristic 
strength of interaction u = U N/K [U [23[ . The quasi- 
linear Rabi regime \u\ < 1, the strong-coupling Josephson 
regime 1 < \u\ < N 2 , and the extremely quantum Fock 
regime u > N 2 . 

With adjustable parameters, the BHH can also real- 
ize an atom interferometer, in which the bias, £, gener- 
ates a phase-shift that can be measured by atom-number 
counting. The interaction term allows for the creation 
of non-classical input states, but also generates unde- 
sired phase-diffusion noise. Atom interferometers based 
on this model are of great current interest because they 
can potentially resolve phase-shifts below the standard 
quantum limit (SQL) of A(f> > 1/^/N, and are limited 
instead by the Heisenberg fundamental limit A(j) > 1/N. 
In such a device, a highly correlated initial state would 



be prepared in the Josephson or Fock regime, but the 
measurement would ideally be made in the Rabi regime. 
While tunable, e.g. via Fcshbach resonance, the interac- 
tion parameter u will never be exactly zero. Thus un- 
derstanding, and potentially harnessing, the dynamical 
effects of the interaction-induced nonlinearity will play a 
crucial role in designing such devices. 

To model the nonlinear dynamics in the Josephson and 
Rabi regimes, it is convenient to define the action- angle 
variables a = e %Vi ^fnl, = y/nie~ lLpi . Using these def- 
initions, the Hamiltonian (TTJ) can be rewritten in terms 
of the relative pair-number n = (rii — 77,2) /2 and relative- 
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In the Josephson interaction regime, states initiated with 
n <C N/2 remain confined to this small population im- 
balance region, so that it is possible to use Josephson's 
approximated Hamiltonian 



H josephson = E c (n - n £ ) 2 - Ej cos(v?) , 



(2) 



with E c = U, Ej = KN/2, and n £ = £/(2U). This 
Hamiltonian matches that of a pendulum with 'mass' 
M = l/(2U) and frequency ujj ps VKUN. The role of 
pendulum angle is played by the relative phase ip between 
the two condensate modes, while the relative population 
imbalance, n, plays the role of angular momentum. One 
consequence of this analogy is that a relative-phase of 
ip = is classically stable (ground state of the pendu- 
lum), whereas a ip — n relative-phase (inverted pendu- 
lum) is classically unstable. 

The term "phase-diffusion" then describes the non- 
linear effects generated by the anharmonic part of the 
cos(p) potential. The many-body manifestation of this 
anharmonicity is the rapid loss of single-particle coher- 
ence for a coherent state prepared with a ir relative- 
phase, as opposed to the slow phase-diffusion of phase- 
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locked condensates with tp> = 



Several recent 



works propose to control such phase-diffusion processes 
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by means of external noise |25|, |26| or modulation of the 
Hamiltonian parameters to induce n flips of the relative 
phase [13. 

In this work, we build further on the pendulum analogy 
to explore the effect of oscillatory driving on the collective 
phase dynamics of the BHH. We consider two possible 
time-dependent driving fields, "vertical" (v) and "hori- 
zontal" (h), given by 



V v , h (t) = f(t)W = D vJl sm(ojt + 



(3) 



Here 'vertical' and 'horizontal' are in reference to the 
pendulum model. The classical phase-pendulum motion 
is in the L X L V equatorial plane of the BHH, with the 
ip = 0, 7r stationary points lying on the L x axis, making 
it the 'gravitation' direction of the pendulum. Hence, V v 
is equivalent to a vertical drive of the pendulum axis and 
Vh corresponds to a horizontal drive. With respect to 
the two- mode BHH, the first type of driving V v (t) is a 
modulation of the hopping frequency K, which may be 
attained by changing the Barrier height, as illustrated 
in Fig. 1(a); whereas Vh(t) may be induced by means 
of shaking the double-well confining potential laterally, 
as illustrated in Fig. 1(b), thus effectively introducing 
the equivalent of an electromotive force in the oscillat- 
ing frame. It is customary to define the dimension- 
less frequency Q = u/uij, and the dimcnsionlcss driving 
strength q = ^UN/K{D/lu) = D/(KCl). Fast and slow 
driving correspond to Q 3> 1 and Q <C 1, respectively, 
whereas q and q <C Q correspond to strong and 

weak driving. 

Within the ID pendulum approximation, the angle 
variable tp and the momentum n are canonical conjugate 
variables. It is well known |28[ that off-resonant, fast 
driving is effectively equivalent to the additional static 
'pseudo-potential ', V° !! h = ±(1 /4)q 2 K£ sin 2 (tp), as illus- 
trated in Fig. 1. It is possible to further refine this effec- 
tive description, adding momentum dependent terms, as 
described in [29[. For sufficiently strong (q 2 > 2) verti- 
cal drive, the effective term V° !! can stabilize the <p = it 
inverted pendulum (Fig. 1(c)), an effect known as the 
Kapitza pendulum [28[. By contrast, the effective term 
V^ ff can destabilize under the same conditions the (p = 
pendulum-down ground state, and generate two new de- 
generate quasi-stationary states (Fig. 1(d)). 

Generally, the driven BHH has a mixed classi cal p hase- 
space structure similar to that of a kicked top [301 ] . with 
chaotic regimes bound by KAM trajectories, making the 
bosonic Josephson junction a good system for studying 
quantum chaos 31, 32| along with the kicked- rotor real- 
ization by cold atoms in periodic optical lattice potentials 
[II] , ultracold atoms in atom-optics billiards [34| , and the 
recent realization of a quantum kicked top by the total 
spin of single 133 Cs atoms [35^ . The unique features of 
the driven BHH in this respect are: (i) It offers a novel 
avenue of 'interaction-induced' chaos, which should be 
distinguished from the single-particle 'potential-induced' 




FIG. 1: (Color online) Schematics of a driven Bose- Josephson 
junction: (a) The Vertical' driving obtained by time- 
dependent modulation of the barrier height between the wells; 
(b) The 'horizontal' driving is via lateral shaking of the 
double-well potential; (c) The potential term in the Joseph- 
son Hamiltonian without driving (dash-dot), and with vertical 
driving (solid) which includes the effective potential (dashed) . 
Circles denote stable stationary points whereas 'x' denotes in- 
stabilities; (d) The same for horizontal driving. It should be 
noted that 'vertical' and 'horizontal' refer to the motion of 
the pendulum axis in the Kapitza analogy, which incidentally 
match the direction of potential modulation 



chaos that had been highlighted in past experiments with 
cold atoms; (ii) The pertinent dynamical variables are 
relative-phase and relative-number, leading to nonlinear 
and possibly chaotic phase dynamics, which may be mon- 
itored via fringe-visibility measurements in interference 
experiments. 

Figure 2 shows representative results for near-resonant 
(Owl) horizontal driving. Stroboscopic Poincare plots 
of the classical (mean-field) evolution at drive-period in- 
tervals are shown on the left for varying drive strength, 
q, demonstrating the growth of the stochastic component 
to form a chaotic 'sea' surrounding regular 'islands' of 
non-chaotic motion. On the right, we represent the full 
many-body BHH evolution via the Husimi Q function 
Q(n,ip) = |(n, tp\ip(t))\ 2 , which provides visualization for 
the expansion of the time-dependent many-particle state 
\if>(i)) in the spin coherent states \n,tp) basis. For weak 
driving the initial preparation |0,0) lies within a regu- 
lar region of phase-space and retains its coherence. In 
contrast, for larger values of q the initial coherent state 
spreads quickly throughout the chaotic sea, resulting in 
a highly correlated many-body state, as manifest in the 

dynamics of fringe visibility (t) = J (L x ) 2 + (L y ) 2 / 1. 
Similar results, with a slightly different classical phase- 
space structure are obtained for vertical driving. 

For an off-resonant weak driving with D, 3> 1 and q <C 
fi, the chaotic component of phase-space becomes hard 
to resolve. In this case we obtain the many-body mani- 
festation of the Kapitza pendulum physics 28[ , where a 
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FIG. 2: (Color online) Mean-field stroboscopic phase-space 
plots of classical trajectories (left), and a representative 
iV=100 many-body quantum Hussimi distribution (right), for 
the dynamics that is generated by the BHH with Vh(t) driv- 
ing, assuming |0, 0) coherent preparation. The parameters 
are it=30, and Q=l. The strength of the driving is q=0.1 
(top), and q=0.5 (middle), and q=1.0 (bottom). The evolu- 
tion of the fringe-visibility is plotted in the lower panel for 
q=0.1 (bold solid, green) and 5=1.0 (solid, blue). The circle 
indicates the time of the Hussimi plot. 



periodic vertical drive of the pendulum axis stabilizes the 
tp = 7r inverted pendulum and a horizontal drive desta- 
bilizes its ground tp = state. 

In generalizing the Kapitza pendulum physics into the 
context of the driven bosonic Josephson junction, we 
have to deal with two modifications: (i) The phase- 
space of the full BHH is spherical and not canoni- 
cal, as opposed to the truncated, cylindrical Joseph- 
son phase-space; (ii) Realistically f(t) may include 
a noisy component whose effect on the effective po- 
tential should be determined. We therefore rc-dcrivc 



the Kapitza physics for the full BHH, using a mas- 
ter equation approach rather than by the standard 
timescale separation methodology. The quantum state 
of the system is represented by the probability matrix 
p, satisfying dp/dt = i[H,p] where W = H+f(t)W with 
f(t) = sm(ojt + <fi). The small parameter is the driving 
period St = 2tt/uj for harmonic drive, or the correlation 
time if f(t) is noisy. Using a standard iterative procedure 
the difference p(t + St) — p(t) can be expressed to 1st or- 
der as an integral over [H+f(t')W,p]. The 2nd order 
adds a double integral over [H+f(t')W, [H+f(t")W,p]], 
and the 3rd order adds a triple integral over 
[H+f(t')W,[H+f(t")W,[H+f(t"')W,p]}]. If /(*) con- 
tains a noisy component, as in the standard master equa- 
tion treatment, we obtain after integration over the sec- 
ond order contribution a diffusion term [W, [W, p]]. In 
the familiar classical Focker-Planck context, with W = x, 
this terms takes the form d 2 p(x , p) / dp 2 . However, for 
a strictly periodic noiseless driving the time integration 
over a period vanishes, and evaluation to 3rd order is 
required. Integrating the 3rd order contribution over a 
period we get terms that can be packed as [[W, [W, H]] , p] . 
Hence, the effective static potential is, 

V'« =~ [W, [W, H]] (4) 

Other terms also exist (producing the tilt of the islands in 
Figs. 3c, 4c), however they depend on the driving phase 
</>, and vanish if the stroboscopic sampling is averaged. In 
the standard canonical case with W = W(x) this expres- 
sion gives the familiar Kapitza result [W (x)] 2 / (AMui 2 ) 
as in (2^ |. For the BHH non-canonical spherical phase- 
space, with W (x L x or W cx L y , it is straightforward to 
verify that Eq.@ generates the expected Kapitza terms 
L 2 or L 2 , as well as additional terms that slightly re- 
normalize the bare values of U and K. 

The predicted Kaptiza physics effects are confirmed 
numerically in Fig. [3] and Fig. 0] for the vertical V v (t) 
and horizontal Vh(t) drive, respectively. Comparison of 
the stroboscopic Poincare plots for the undriven (a) and 
driven (c) BHH, clearly shows the stabilization of the 
|0,7r) coherent state by the vertical L x driving (Fig. [3]), 
and the destabilization of the |0, 0) preparation by the 
horizontal L y driving (Fig. [4]). These classical effect are 
mirrored in the evolution of the quantum Husimi func- 
tion, thereby affecting the many-body fringe-visibility 
dynamics, leading to the protection of coherence by V v {t) 
driving for tp — ir coherent preparation, and to its de- 
struction by Vh(t) for tp = 0. 

To conclude, the driven BHH, currently attainable in a 
number of experimental setups, presents a wealth of non- 
linear phase-dynamics effects. Strong, resonant driving 
fields result in large chaotic phase-space regions, open- 
ing the way for t he g eneration of exotic highly correlated 
quantum states [3l|. The properties of such states, as 
well as their manifestation in interference experiments, 



4 



50 



-50 
50 



-50 



Jl 



2ji 





FIG. 3: (Color online) Quantum Kapitza pendulum with ver- 
tical driving, and |0, n) coherent preparation. The panels are 
arranged as in Fig. 1. The parameters are u=100, fi=30, and 
iV=100. The strength of the driving is q=0 (a,b, lower panel 
dashed blue) and q—S (c,d, lower panel solid red). Circles 
denote the time at which the Husimi distribution in (c,d) is 
plotted. The classical stabilization of the inverted pendulum 
results in a protected single-particle coherence of the initial 
preparation. 



and the more conventional tunneling effects between reg- 
ular islands, are novel manifestations of semiclassical 
physics. For weak and fast off-resonant drive we have ob- 
tained the many-body equivalents of the Kapitza pendu- 
lum effects, with the relative-phase between the conden- 
sates acting as the pendulum angle. Such effects could be 
readily observed in interference experiments and utilized 
to protect fringe-visibility. We note that noise-protected 
coherence was also studied in Ref . j2f| , yet with a rather 
different quantum-Zcno underlying physics. 
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FIG. 4: (Color online) Quantum Kapitza pendulum with hor- 
izontal driving, and |0,0) coherent preparation. The panels 
are arranged as in Fig. 3, with the same parameters, except 
u = 30. One observes the splitting of the ground state. Due 
to this destabilization the fringe visibility of the preparation 
is destroyed. 
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